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NMR SPECTROSCOPY USING SPIN STATES WITH LONG LIFETIMES 

The present invention relates to the field of NMR. More specifically, the present 
invention relates to a method of NMR spectroscopy which utilises non-equilibrium spin 
states with long lifetimes. 

Nuclear magnetic resonance (NMR) is a powerfiil and well established tool for studying 
samples and sample interactions. It can also be used for imaging, for example. 
Magnetic Resonance Imaging (MRI). 

In NMR, the spin and magnetism of atomic nuclei are exploited to provide inforaiation 
about the chemical composition, spatial distribution, or molecxilar motion of the 
molecules or atoms. 

One of the limitations of NMR is the low intrinsic signal strength. Attempts to address 
this problem have included nuclear spin hyperpolarization experiments, for example 
using parahydrogen (C. R. Bowers and D. P. Weitekamp, Phys. Rev. Lett. 57, 2645- 
2648 (1986); C. R. Bowers and D. P. Weitekamp, J. Am. Chem. Soc. 109, 5541-5542 
(1987); M. G. Pravica and D. P. Weitekamp, Chem. Phys. Lett. 145, 255-258 (1988)). 

One of the principal factors contributing to the applicability of NMR based techniques 
is the relatively long lifetime of the nuclear spin order. This long lifetime allows NMR 
to be used to follow processes such as diffusion, flow, slow molecular motion and 
chemical reactions. The relaxation of the nuclear spins back to thermal equilibrium, and 
hence the lifetime of the nuclear spin memory is characterised by a time constant, Ti 
known as the longitudinal relaxation time constant or as the spin lattice relaxation time 
constant. The time constant Ti is governed by molecular motion, and is in general 
different fi-om the phase memory time T2, which depends upon molecular motion and 
also external perturbations such as magnetic field gradients. In general, Ti is always 
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longer than O.5T2, although in common situations, Ti is much greater than T2, often by 
an order of magnitude. 

Generally, for many samples, the time constant Ti is of the order of seconds. Thus, the 
time for following a process of the type meationed above is very limited. 

Attempts have been made to extend the phase memory tune in the context of quantum 

computation, using the concept of Decoherence-Free Subspaces (L. Viola, E. M. 

Fortunato, M. A. Pravia, E. Knill, R. Laflamme and D. G. Cory, Science 293, 2059- 
2063 (2001); E. M. Fortunato, L. Viola, M. A. Pravia, E. Ml, R. Laflamme, T. F, 
Havel and D. G. Cory, quant-ph /0210057 (2002); D. A. Lidar and K. B. Whaley, quant- 
ph /0301032 (2003); J. E. OUerenshaw, D. A. Lidar and L. E. Kay, quant-ph /0302175 
(2003)). However, although the compensation of artificial noise contributions to the 
phase memory time T2 has been demonstrated, the extension of the spin memory time 
beyond Ti has not been demonstrated. 

The present inventors have found that it is possible to generate nuclear spin states which 
remain far from equilibrium for a time much longer than Ti. These spin states may be 
advantageously used to study processes witii a longer lifetime than that of Tj. These 
non-equilibrium spin states may be studied (and in some cases, created) by switching 
the symmetry of the nuclear spin Hamiltonian operator. 

Thus, in a first aspect, the present invention provides an NMR method, the method 

comprising: 

providing a sample where the nuclear spin Hamiltonian operator of the component 

molecules of the sample possess one or more symmetry operations; 

creating a quasi equilibrium nuclear spin ensemble state in a sample, said quasi 

equilibrixmi nuclear spin ensemble state comprising at least two manifolds of spin states 

which transform differently under said symmetry operations of said Hamiltonian and 

said manifolds having different mean nuclear spin populations; 

allowing said quasi equilibrium nuclear spin ensemble state to remain for a time of 

equal to or substantially greater than 3Ti, where Ti is the spin lattice relaxation time; 
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breaking at least one of the said symmetry operations of said Hamiltonian; 

applying a sequence of magnetic fields to generate a nuclear magnetic resonance signal 

from said sample; and 

detecting said nuclear magnetic resonance signal. 

The temi symmetry is used here to describe an exchange operation of one or more of the 
nuclear spins that leaves the nuclear spin Hamiltonian invariant. In practice, ttie 
imposition of a new symmetry operation does not need to be exact. This is discussed in 
more detail later. It will be appreciated by those skilled in the art that the said 
Hamiltonian is the effective Hamiltonian of the nuclear spins, which governs the 
evolution of the nuclear magnetic moments over the relevant part of the appUed 
experimental sequence. 

Manifold is the term used in the art for groups of nuclear spin states that transform in a 
common way under the symmetry operations of the nuclear spin Hamiltonian, The 
manifolds may reach internal equilibrium independently and may have different mean 
nuclear spm populations. A state of intemal equiUbrium for the individual manifolds 
but different mean nuclear spin populations for different manifolds is, by definition, not 
a thermal equilibrium state for the entire ensemble of nuclear spins. Nevertheless, it 
may be very long lived, since there are no relaxation processes mternal to the molecules 
which can lead to its destruction, hence it is referred to as a quasi equilibrium nuclear 
spin ens^ble state. 

For the avoidance of doubt and as used herein, the term quasi equilibrium nuclear spin 
ensemble state will be used to mean a state of the spin ensemble in which manifolds of 
different symmetry have reached intemal equilibrium but in which the different 
manifolds are not equilibrated with each other. 

An example of a quasi-equiUbrium spia state is foxmd in an ensemble of molecules each 
containing two spins-1/2 of the same isotopic type, in a magnetic field which is 
sufficiently low that the chemical shift differences between the nuclei may be ignored. 
In this case, an exchange of the two nuclear spms is as a symmetry operation of the 
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nuclear spin Hamiltonian operator. The nuclear spin states may be classified as 
symmetric with respect to this operation (the three triplet states) or as antisymmetric 
with respect to fliis operation (the singlet state). A state of the ensemble of molecules, 
in which the populations of the triplet states have equilibrated wifli each other, but in 
which the singlet population is different fiom the mean triplet population, qualifies as a 
"quasi-equilibrium nuclear spm state". 

Another example is found in a sample composed of molecules of the form A-CRt- . 
CHr-X, where A and X are arbitrary groups, supposed to have no significant nuclear 
spin magnetism. In low magnetic field, an exchange of the spins in the two CH2 groups 
with respect to each other is a symmetry operation of the nuclear spin Hamiltonian 
operator. The nuclear spin states may be classified as symmetric, or antisymmetric, with 
respect to this spin exchange operation. A state of the ensemble of molecules, in which 
the populations of the symmetric states have equilibrated with each other, and in which 
the populations of the antisymmetric states have eqiulibrated with each other, but in 
which the mean symmetric state population is different firom the mean antisymmetric 
state population, qualifies as a "quasi-^equilibrium nuclear spin ensemble state*'. 

A non-equilibrium singlet population is a special case of a quasi-equilibrium spin state, 
which applies when the molecules ovly contain two nuclei with significant magnetic 
properties. Quasi-equilibrium nuclear spin ensemble states apply to molecules with two 
or more coupled nuclear spins. The present invention may be used witii a sample which 
can form a non-equilibrium singlet population or with a sample which can form a quasi- 
equilibriiun spin ensemble state other than a non-equilibrium singlet population. 

In a second aspect, the present invention provides an NMR method, the method 
comprising: providing a sample where the nuclear spin Hamiltonian operator of the 
component molecules of the sample possess one or more symmetry operations; creating 
a quasi equilibrium nuclear spin ensemble state in a sample, said quasi equilibrium 
nuclear spin ensemble state comprising at least two manifolds of spin states which 
transform differently xmder said symmetry operations of said Hamiltonian and said 
manifolds having different mean nuclear spin populations, said quasi equiUbrium 
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nuclear spin ensemble state not having a singlet state population; breaking the 
synnnetry operation of said Hamiltonian; applying a sequence of magnetic fields to 
generate a nuclear magnetic resonance signal fix>m said sample; and detecting said 
nuclear magnetic resonance signal. 

In the method of the second aspect of the present invention, the quasi-equiUbrium 
nuclear spin ensemble state lasts for a time of at least 3Ti, Thus, the method may 
further comprises allowing said quasi. equiUbriuin nuclear spiii ensemble state to remain 
for a time of equal to or greater than 3Ti, where Ti is the spin lattice relaxation time 
before breaking at least one symmetry operation of said Hamiltonian; 
In any of the above aspects of the present invention, a sample is prepared which has a 
Hamiltonian possessing one or more symmetry operations, this state can exist for a time 
of at least 3Ti before one or more of the symmetry operations is removed from the 
Hanoiltonian. 

A number of different methods may be used to impose and remove a symmetry 
operation for the nuclear spin Hamiltonian operator and some of these are dependent on 
the nature of the sample. 

When the molecules of the sample comprise nuclear spins of different isotopic types or 
have different chemical shifts and the molecules have sufficient geometrical symmetry 
said step of providing a sample where the nuclear spin Hamiltonian operator of the 
component molecules of the sample possess one or more symmetry operations may 
comprise reducmg a magnetic field to a value where a symmetry operation is imposed 
on said Hamiltonian. A molecule of sufficient geometrical symmetry is defined as a 
molecule for which the spin Hamiltonian, including all direct and indirect spin-spin 
couplings, has higher symmetry when the chemical shifts are ignored than when the 
chemical shifts are present. 

Altematively, the sample could be kept in hi^ field while a radiofrequency pulse 
sequence is appUed which has the effect of suppressing the chemical shift differences. 
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see for example J. D. EUett and J. S. Waugh, J. Chem. Phys. 51, 2851 (1969); M. H. 
Levitt, R. Freeman and T. A. Frenkiel, Adv. Magn. Reson. 11, 47 (1983). 

Alternatively, a sample may be subjected to a chemical reaction, configured to switch 
the molecules between different conformational or configurational states, in which spin 
symmetry operations are imposed in order to provide a sample where the nuclear spin 
Hamiltonian operator of the component molecules of the sanaple possess one or more 
symmetry operations. An example of such a chemical reaction is a photoisomerization 
reaction. Such a reaction may be used to change a molecule firom a mixed cis-trans 
form to a pure cis or pure trans form. 

A number of different methods may be employed for creating the quasi-equilibrium 
nuclear spin ensemble state. For example, the sample maybe subjected to 
radiofrequency pulses of magnetic field in a high magnetic field environment. 
Specifically, thermal equilibration of a sample in high magnetic field, perturbation of 
the eqmlibrium of the nuclear spm system using apphed radiofirequency pulses, and 
transport of the sample into the low magnetic field. 

Those skilled in the art are aware of the technique of manipulating the populations of 
spin states using timed sequences of pulsed radiofirequency fields, resonant with the 
precession fi:equency of the nuclei in the high static magnetic field. The sequence of 
pulsed fields may be configured to produce a non-equilbrium state in high field, which 
gives rise to a quasi-equiUbrium configuration of spin state populations, once the 
sample is transported to low magnetic field. 

The quasi-equilibrium nuclear spin ensemble state may also be provided by reacting an 
unsaturated symmetric molecule with parahydrogen. 

Alternatively, the quasi-equilibrium nuclear spin ensemble state may also be provided 
by reacting an unsaturated non-symmetric molecule with parahydrogen, the product of 
said reaction being placed in a magnetic field which is sufficiently low such that the 
difference in chemical shifts be ignored. 
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Providing a sample with the desired symmetry and creating a quasi-equilibrium nuclear 
spin ensemble state, may be performed as distinct steps in sequence or otherwise, or 
may be performed together, for example by manipulating the populations of the spin 
states of a molecule such that a quasi-equilibrium spin state is formed once a symmetry 
operation is provided to the Hamiltonian. 

Since there are no relaxation processes internal to the molecule which can convert the 
quasi-equiUbrium nuclear spin ensemble state to a fully equiUbrated ensemble state, the 
lifetime of the quasi-eqmlibrimn nuclear spin ensemble state can be much longer than 
the conventional relaxation time Ti.. For this reason, such states may be used a long- 
term repository of nuclear spin order. For example, the inventors found that, in systems 
of inequivalent spin pairs in low magnetic field, the decay of the singlet state 
population, or other quasi-equiUbrium nuclear spin ensemble state, is governed by a 
time constant which is much longer than the spin relaxation time Ti This long lifetime 
allows nuclear spin order to persist during processes which take longer than Ti. 

Thus, while said quasi-equilibrium nuclear spin ensemble state is maintained, the 
sample may flow through an apparatus or a body, may be reacted with a second sample, 
may be purified or placed in contact with a second sample, in order to allow a transfer 
of nuclear spin order. The quasi-equilibrium nuclear spin ensemble state is preferably 
maintained for at least 5Ti, more preferably at least lOTi. 

The quasi-equilibrium nuclear spin ensemble state may be observed by (i) breaking the 
symmetry operation of the molecular spin Hamiltonian followed by (ii) a sequence of 
events leading to the genoation of an informative NMR signal. 

The operations leading to the breaking of symmetry are the same as those discussed 
above in the context of the making of symmetry. For example, in systems displaying 
chemical ineqmvalence, the symmetry may be broken by imposing a sufficiently large 
external magnetic field. In other cases, the symmetry may be broken by 
photoisomerization processes or other chemical reactions. 
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Applying a sequence of magaetic fields to generate an NMR signal from said sample 
may be perfonned in eith^ a high static magnetic field backgromid or a low static 
magnetic field background. A low field is considered to be a field which is equal to or 
less than 20mT. 

The condition for a low field is that the difference in the precession firequencies of the 
involved spins must be equal to, or snialler than the J^couplings. In the case of two 
protons, with a chemical shift difference of lOppm and a J-coupling of 1 Hz (roughly 
worst case), this condition implies that the field should be less than around 2mT, which 
is easy to achieve, even without magnetic shielding. If the chemical shift difference is 
smaller, or if the J-coupling is larger, larger magnetic fields may be tolerated, in direct 
proportion to the J-coupling, and in inverse proportion to the chemical shift difference. 
The heteronuclear case is more difiBcult; for example, for a system of protons and 
carbon-13 nuclei, with a J-coupling of 120 Hz, the external field must be less than 3 jiT, 
which requires careftil magnetic shielding. A high field condition as used herein is 
when the difference in the precession firequencies of the involved spins is larger than the 
J-couplings. 

The NMR signal may be detected in the standard way using electromagnetic induction 
within a coil. Alternatively, a direct magnetic field detecting device such as 
Superconducting Quantum Inteference Device "SQUID" may be used (see for example 
R. McDermott, A. H. Trabesinger, M. Mtlck, E. L. Hahn, A. Pines and J. Clarke, 
Science 295, 2247-2249 (2002)). 

In a third aspect, the present invention provides an NMR apparatus, the apparatus 
comprising: means for providing a sample where the nuclear spin Hamiltonian operator 
of the component molecules of the sample possess one or more symmetry operations; 
means for creating a quasi equilibrium nuclear spin ensemble state in a sample, said 
quasi equilibrium nuclear spin ensemble state comprising at least two manifolds of spin 
states which transform differently under said symmetry operations of said Hamiltonian 
and said manifolds having different mean nuclear spin populations, said quasi 
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equilibrium nuclear spin ensemble state not having a singlet state population; means for 
breaking the symmetry operation of said Hamiltonian; means for applying a sequence of 
magnetic fields to generate a nuclear magaetic resonance signal from said sample; and 
means for detecting said nuclear magnetic resonance signal. 

The means for providing a sample where the nuclear spin Hamiltonian operator of the 
component molecules of the sample possess one or more symmetry operations may be 
provided by a chemical or physical enyiromnent coiifigured to produce the desired 
sample. For example, a means may be provided to perform any of the above suggested 
methods for providing such a sample. 

Similarly, means for creating a a quasi equilibrium nuclear spin ensemble state in said 
sample may be a means for allowing performance of any of the above described 
methods for creating the a quasi equilibrium nuclear spin ensemble state. 

Although the above discussion has concentrated on an NMR method, the creation of a 
quasi equilibrium nuclear spin ensemble state may be used for other applications, for 
example for storing information in the form of quantum bits or '*qubits". 

Thus, in a fourth aspect, the present invention provides a method of storing information, 
said method comprising: providing a sample where the nuclear spin Hamiltonian 
operator of the component molecules of the sample possess one or more symmetry 
operations; and creating a quasi equihbrium nuclear sphi ensemble state in a sample, 
said quasi equihbrium nuclear spin ensemble state comprising at least two manifolds of 
spin states which transform differently under said symmetry operations of said 
Hamiltonian and said manifolds having different mean nuclear spin populations; 
assigning a bit value to at least one spin state of said quasi equihbrium nuclear spin 
ensemble state and storing said spin state for a time equal to or substantially greater than 
3Ti, where Ti is the spin lattice relaxation time. 

In a fifth aspect, the present invention provides a method of storing information, said 
method comprising: providing a sample where the nuclear spin Hamiltonian operator of 
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the componeat molecules of the sample possess one or more symmetry operations; and 
creating a. quasi equilibrium nuclear spin ensemble state in a sample, said quasi 
equilibrium nuclear spin ensemble state comprising at least two manifolds of spin states 
which transform differently under said symmetry operations of said Hamiltonian and 
said manifolds having different mean nuclear spin populations, said quasi equilibrium 
nuclear spin ensemble state not having a singlet state population; assigning a bit value 
to at least one spin state of said quasi equilibrium nuclear spin ensemble state and 
storing said spin state for a time equal to or substantially greater than 3Ti, where Ti is 
the spin lattice relaxation time. 

The preferred methods discussed above in relation to NMR for providing the sample 
and creating the quasi equiUbrium nuclear spin ensemble state may also be used with 
the above fourth and fifth aspects of the invention. 

The above methods may be used to store qubits in quantum computing or related 
applications. Previously, NMR techniques have been suggested for quantum 
computation by: D. G. Cory, A. R Fahmy and T. F. Havel, Proc. Natl. Acad. Sci. USA 
94, 1634-1639 (1997); N, A. Gershenfeld and 1. L. Chuang, Science 275, 350-356 
(1 997). The techniques covered in these papers may be combined with the present 
invention. 

In a sbcth aspect, the present invention provides an ^paratus for storing information, 
said apparatus comprising: means for providing a sample where the nuclear spin 
Hamiltonian operator of the component molecules of the sample possess one or more 
symmetry operations; and means for creating a quasi equilibrium nuclear spin ensemble 
state in a sample, said quasi equiUbrium nuclear spin ensemble state comprising at least 
two manifolds of spin states which transform differently under said symmetry 
operations of said Haimltonian and said manifolds having different mean nuclear spin 
populations; means for assigning a bit value to at least one spin state of said quasi 
equiUbriimi nuclear spia ensemble state and storing said spin state for a time equal to or 
substantially greater than 3Ti, where Ti is the spin lattice relaxation time. 
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In a seventh aspect, the present invention provides an apparatus for storing information, 
said apparatus comprising: means for providing a sample where the nuclear spin 
Hamiltonian operator of the component molecules of the sample possess one or more 
synunetry operations; and means for creating a quasi equilibrium nuclear spin ensemble 
state in a sample, said quasi equilibrium nuclear spin ensemble state conq)rising at least 
two manifolds of spin states which transform differently under said symmetry 
operations of said Hamiltonian and said manifolds having different mean nuclear spin 

populations, said quasi equilibrium nuclear spin ensemble state not haying a singlet 

state population; means for assigning a bit value to at least one spin state of said quasi 
equilibrium nuclear spin ensemble state. 

The present invention will now be described with reference to the following preferred 
non-limiting embodiments in which: 

Figure 1 is a flow diagram outlining a method in accordance with an embodiment of the 
present invention; 

Figure 2 is a flow diagram to illustrate the switching of symmetry operations in the 
Hamiltonian; 

Figure 3 is a flow diagram illustrating a method in accordance with an embodunent of 
tiie present invention; 

Figure 4a is a schematic plot of B-field applied to the sample against time, Figure 4b is 
a schematic plot of radiofirequency magnetic pulse sequences applied perpendicular to 
the sample against time, and Figure 4c is a schematic plot of the distribution of spin 
state populations of a sample subjected to the static fields of Figure 4a and the pulse 
sequences of Figure 4b; 

Figure 5a is a conventional ^H NMR spectrum of the solution of 2,3-dibromothiophene 
(inset) in DMSO-d^, Figure 5b is a spectrum generated by the method of Figures 3 and 
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4a, and Figure 3c is a spectrum generated by the method of Figures 3 and 4a, using a 
diflferent pulse sequence to t used to produce the spectrum of Figure 5b; 

Figure 6 schematically illustrates a plot of signal amplitude against time when the 
sample is placed in a low field condition; 

Figure 7 are plots of singlet population against time for a number of different magnetic 
-fields-which are to be used as the low field,.Figure 7a is for a lo\y field of OT, Figure 7b 
is for a low field 2QmT, Figure 7c is for a low field of 200mT, Figure 7d is for a low 
field of 50QmT, Figure 7e is for a low field of IT and Figure 7f is for a field of 5T; and 

Figure 8 is a flow diagram schematically illustrating a method in accordance with a 
fiirfher embodiment of the present invention; 

Figure 9 is a diagram schematically illustrating a chemical reaction of a symmetric 
molecule with parahydrogen; 

Figure 10 is a flow diagram schematically illustrating a method in accordance with a 
fiirther embodiment of the present invention; 

Figure 1 1 is a diagram schematically illustrating a chemical reaction of a non-symmetric 
molecule with parahydrogen; 

Figure 12 is a flow diagram schematically illustrating a method in accordance with a 
fijrfher embodiment of the present invention; 

Figures 13a and 13b schematically illustrate a photoisomorization process which may be 
used to impose or remove a symmetry operation on a Hamiltonian; and 

Figure 14 is a flow diagram schematically illustrating a method in accordance with a 
further embodiment of the present invention. 
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Figure 1 outlines the basic steps in accordance with an embodiment of the present 
invention. 

In step SlOl, a sample is provided where the nuclear spin Hamiltonian operator of the 
component molecules of the sample possess one or more symmetry operations. This 
step may involve providing a sample which has this quality inherently (such a method 
will be described in more detail with reference to figure 8), alternatively, the sample 
may only possess this quality in certain environments (this \vill be described with 
reference to figures 3, 10 and 14), as a fiirther option, the sample may need to be 
chemically modified to achieve the symmetry requiremrat (this will be described with 
reference to figure 13). 

In Step S103, a quasi-equilibrium nuclear spin ensemble state is created, how this state 
may be created will be described later with reference to figures 3 to 14. The term quasi 
equilibrium nuclear spin ensemble state is used to refer to the state of a spin ensemble in 
which manifolds of different symmetry have reached intemal equilibrium but in which 
the different manifolds are not equilibrated with each other. 

Once the quasi equilibrium nuclear spin ensemble state has been created, it may be 
stored for a reasonably long time, much longer the spin-lattice relaxation time Ti . Ti is 
believed to be the time constant which determines when the spin memory is quenched 
in NMR experiments. However, the inventors have found that the quenching of quasi 
equilibrium nuclear spin ensemble states is govemed by a time constant Tqe which is 
much greater than Ti. This allows the sample to undergo processes SI 05 before the 
quasi equilibrium nuclear spin ensemble state decays and for a time greater than or 
equal to 3Ti. For example, the sample may be injected or inserted into a human or 
animal body and its passage through the body may then be studied vising NMR. Since 
the sample may be allowed to propagate in the body for longer than conventionally 
allowed due to the discovery of the long lifetime of the quasi equilibriimi nuclear spin 
ensemble state, more information may be gathered and different types of experiments 
may be conducted utiUsing the long lifetime of the quasi equilibrium nuclear spin 
ensemble state. 
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The extension of spin memory lifetimes beyond Ti by using quasi-equilibrimn states is 
distinct from the extension of spin dephasing lifetimes beyond T2 demonstrated in the 
context of Decoherence-Free Subspaces as employed in quantum computation (L. 
Viola, E. M. Fortunato, M. A. Pravia, E. Knill, R. Laflamme and D. G. Cory, Science 
293, 2059-2063 (2001); E. M. Fortunato, L. Viola, M, A. Pravia, E. Knill, R, Laflamme, 
T. F. Havel and D. G. Cory, quant-ph /0210057 (2002); D. A. Lidar and K. B. Whaley, 
quant-ph /0301032 (2003); J. E. Ollerenshaw, D. A. Lidar and L. E. Kay, quant-ph 
• /0302175 (2003))..The main distmctions are: (i) the major Ti mechanism is usually the 
modulation of the dipole-dipole coupling by random molecular motion, while the 
experiments reported in the quantum computation field demonstrate the compensation 
of an artificial shortening of T2 by molecular diffiision in strong magnetic field 
gradients; (ii) the methods reported in the above references are incapable of 
demonstrating a lengthening of the spin memory time beyond Ti. 

The sample may be reacted with other samples or may undergo a purification process 
during step SI 05. Alternatively, the sample may be placed in physical contact with 
other molecules to transfer polarisation of the quasi equilibrium nuclear spin ensemble 
state. 

The quasi equilibrium nuclear spin ensemble state may only be observed if the 
symmetry of the sample is broken S107. This may be done by a number of methods 
which will be described with reference to figures 2 to 13. 

The sample is then subjected to a sequence of magnetic fields to generate nuclear spin 
magnetization. The sequence of magnetic fields may be ^plied in a high magnetic 
field environment as described in detail with reference to Figure 4, or a low or zero 
magnetic field envuromnent as described in R. McDermott, A. H. Trabesmger, M. 
Miick, E. L. Hahn, A. Pines and J, Clarke, Science 295, 2247-2249 (2002). 

The signal from the NMR active sample can then be detected using either a high field 
detection technique such as electromagnetic induction of current in a coil or direct 
detection of the magnetisation using a SQUID or the like. 
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Figure 2 schematically illustrates the symmetry considerations for steps S103 to S109 of 
figure 1. The small rectangles represent two different state manifolds, and the shading 
represents their populations. In S201, the quasi equilibrium nuclear spin ensemble state 
has been prepared and the sample has a Hamiltonian which has one or more symmetry 
operations. In this enviroimient, the quasi equilibrium nuclear spin ensemble state 
decays with a long time constant Tqe, which is much longer than the spin lattice 
relaxation tune Ti. . 

In step 203, the symmetry operations are removed from the Hamiltonian, allowing the 
NMR signal to be observed,. Due to spin-lattice relaxation, the spin population between 
the states equilibrates and an equilibrium state is reached after Ti in step S205. 

Figure 3 and figure 4 schematically illustrate a method in accordance with an 
embodiment of the present invention. Figure 3 is a flow diagram summarising the 
processes, while Figures 4a, 4b and 4c describe the changes to the sample and its 
environment. 

In NMR procedures, a sample is subjected to a high constant magnetic field Bhigh* A 
much smaller rf pulsed magnetic field is then provided perpendicular to the high 
magnetic field Bhi^. 

The pulsed magnetic field follows a puke sequence A+. The pulse sequence is applied 
at the Larmor firequency of the spins in the high magnetic field. The pulse sequence 
will be described in greater detail later in the description. It is well known to those 
skilled in the art of NMR how to design pulse sequences in order to achieve a desired 
distribution of spin states within a sample. For the purposes of symmetry-breaking, the 
high magnetic field simply needs to be an order of magnitude larger than the "low" field 
discussed above, i.e. more tfian a few mT for typical proton systems, and more than a 
few nT for heteronuclear systems. However, for the purposes of sensitive observation of 
the NMR signals, the field will usually be as large as possible, i.e. between 1 and 20 
Tesla, depending upon the instrument available. 
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Figure 4c schematically illustrates how the spin states may be manipulated using pulse 
sequence Af . 

At time point 1, the nuclear spin system has four Zeeman energy levels with energies in 
the order Eaa (3) < Eap (5) < Epa (7) < Epp (9). In the figure, Eap (5) and Epa (7) are 
shown to be virtually identical. However, there is a slight energy difference between 

these two levels. The above assumes that the chemical shift .82 is greater than the 

chemical shift 81, the gyromagnetic ratio y and the J-coupling are both positive and that 
the spin system is weakly-coupled in high fields. The population is distributed 
throughout the foxir states. In thermal equilibrium at ambient temperature, the lower 
most level | aa> has a small excess population (depicted by the filled balls), while the 
uppemiost level | pp> has a slightly depleted population (depicted by the white balls). 
These population differences are very small in the case of conventional fliermal 
polarisation, but may be of the order of one in spin hyperpolarization experiments. The 
sample is in thermal equilibrium and this stage of the sample corresponds to step S301 
in figure 4. 

The pulse sequence A+ is then appUed until time point 2, taking the sample to step S303 
in figure 4. During this time, the population of state | aa> exchanges with that of 
I aP>, and that of | poO with that of | PP>. Thus, a non-equiUbrium population 
distribution is generated, in which level 1 ap> has an excess population, while | Poc> 
has a depleted population. This is shown at time point 2. 

During the interval from time point 2 to time point 3, the sample is transported to 
a region of low magnetic field Biow, which is of the order of milliTesla, taking the 
sample to S305 of figure 4. Typically, Blow wiU be at most 20mT. hi low field, it is 
convenient to discuss the behaviour of the nuclear spins using the four eigenstates of the 
pure J-coupling Hamiltonian [1]. 

Hj^lnJI,^!^ [1] 
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The four J-cotqpling eigenstates may be classified as the three components of a triplet 
state, plus a smglet state. The triplet state components are denoted by equation [2]. 

l?;)=l««> 

\T,)=^i\ap)^\Pa)) [2] 



The triplet states have the same J-coupling eigenvalue of +-1 tiJ. The singlet state is 
denoted by equation [3]. 

\S,)=j=i\aP)-\pa)) [3] 

The singlet state has a J-coiq)ling eigenvalue of -f ;c/ , The triplet and singlet states 

are energy eigenstates of the nuclear spins only if the magnetic field is zero. The energy 
level diagram in this case is depicted next to time point 3 m Figure 4c. Although v/e use 
the zero-field eigenstates as a basis for a representation of the spin density operator, the 
treatment below does not assume a field of exactly zero. 

If the transport from high to low field is fast compared to Ti but slow compared to the J- 
coupling, the excess population of the high-field eigenstate | ap> is transferred 
adiabatically into an excess population of the low-field eigenstate | So>. Similarly, the 
depleted population of the high-field eigenstate | Pa> is transferred adiabatically into a 
depleted population ofthe zero-field eigenstate |To>. The idealized distribution of the 
singlet and triplet state populations is shown at time point 3 in Figure 4c. This assumes 
pure adiabatic transfer and neglects Ti relaxation during the transport interval. 

The nuclear spins evolve in low magnetic field under a combination of coherent and 
incoherent relaxation processes. Here we use a standard second-order perturbation 
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treatment of the modulate dipole-dipole interaction to analyze the spin-lattice relaxation 
in low field. 

In general, the spin density operator p evolves according to the Liouville-von Neumann 
equation [4]. 

^p{thLp(t) [4] 

Where the low-field Liouvillian is given by equation [5]. 

4=tr-i^^ [5] 

H^j, is the superoperator for the commutator with the low-field spin Hamiltonian of 
equation [6], 

Hjj. = -rB,ow (l + ^1 Viz ~ /Blow (l + ^2 + Hj [6] 

The thermal polarization terms may be ignored in low field. If the relaxation is 
domiiiated by the intramolecular dipole-dipole coupling mechanism, the relaxation 
superoperator may be evaluated by second-order time-dependent perturbation theory 
and is given by equation [7]. 



f^=-£i^^(-ir)^^^(0Vr [7] 

^DD (0 is ^® commutation superoperator for the dipole-dipole coupling Hamiltonian, 
evaluated at time t, and the overbar represents an ensemble average. Equation [4] is 
general for any field strength low enough to satisfy tiiie extreme narrowing condition for 
the dipole-dipole relaxation. 



18 



wo 2005/015253 



PCT/GB2004/002966 



It is convenient to ns the matrix Tepiesentations of the Liouvillian siq>eToperator Z^^^ in 
a space of orthonoimal spin operators. The following set of six operators generate a 
suitable basis: 

|1) = ISoxSol =i(/r/f+/f/?-ir/2-/,-i^) 

|2) = iTiXTil =1^1^ 

|3) ^ iToXTol =^{i^il+ifit+i:i;+M.J:) . . 

|4) = |T-i><T.,l =//'/f 

|5) =i(i5oXrol+|ro){5o|) =^(A"/f -itn) 



where /t are single-spin shift operators and {/ J /are single-spin polarization 
operators. Theket 1 1) describes the population operators of the singlet state. Kets U), 
1 3) and 1 4) describe the population operators of the three triplet states. Kets 1 5) and 
1 6) represent the x and y-components of the coherence between the singlet state and the 
cmtral component of the triplet. All of these operators commute with the angular 
momentum operator along tiie field axis, as does the low-field Hamiltonian Hjj? . The 
operator basis of equation [8] is closed under the appUcation of the low-field Liouvillian 

The matrix representation of the low-field relaxation superoperator f^^ may be 
evaluated by standard relaxation theory, assuming pure intramolecular dipole-dipole 
relaxation driven by rigid isotropic random rotation with a correlation time Xc. The low- 
field Liouvillian evaluates to 
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[8] 

where the dipole-dipole coupling hetween the spins is given by b= -(|io/47c)y^fir'^, and r 
is the intemuclear distance. The low-field chemical shift firequency diflFerence is 
defined by equation [9] . 

In zero field o)^^ = 0 , the first row and column of equation [8] contains only zeros. 
This indicates that the singlet state population is dynamically isolated in low field, and 
is conserved under the evolution process. The zero-field singlet state may be used as a 
long-term repository for nuclear spin order, protected fix)m intramolecular dipole-dipole 
relaxation. 

If tihie magnetic field is not exactly zero, equation [8] indicates that the singlet state is 
relaxed indirectly by dynamic couplings to the central triplet state and to singlet-triplet 
coherence. 

The idealized behaviour of the spin ensemble after an interval tlf in low field is 
depicted by the diagram at time point 4 in Figure 4c. If the field is sufficiently low, and 
Tlf is long enough, the triplet populations equilibrate, while the singlet population is 
conserved. Thus, a quasi-eqmUbrium ensemble spin state is rapidly established. This 
quasi-equilibrium ensemble spin state is stable for a time much longer than Ti. 

The singlet population may be read out be adiabatic transport of the sample back into 
high field during the interval z*®^. This process breaks the symmetry of the sample 
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and forms a sample which is an active NMR sample as illustrated at step S307 of figure 
3. 



Specifically, this process transforms the low-field singlet population into a population 
of the high-field | aP> state, as depicted for time point 5 in Figure 4c. The population 
of the I ap> state is conv^ted into observable NMR signals by the pulse sequence B at 
step S309 of figure 3. 

Sequence B is given by three rf pulses separated by two unequal delays, and converts 
the 1 ap> state population into antiphase NMR signals centered at the chemical shift 82. 
NMR signals deriving fi-om Ti relaxation during tiie transport interval, on the other 
hand, only generate NMR signals at the chemical shift 5i. The interesting NMR signals 
derived from low-field singlet storage are therefore cleanly separated firom signals with 
a trivial origin. 



The pulse sequence A+ is given by the two-pulse sequence 90o - xi - 9O.90. where the 
symbol p^, denotes a strong, non-selective pulse with flip angle P and phase <|) (both 
angles are specified in degrees, and the phases taken into account the sign of the 
precession and the radio-firequency mixing scheme. The delay xi is set to the value 




where 



<"'=-)^U^2-sd [10] 



For the experimental results, the delay was set to the value xi = 2.05ms. 

Standard spin operator theory shows that in the case 62 > 5i and y > 0, the pulse 
sequence Af has the ^proximate propagator 



C/Oi)= exp[-z^(/^+/2,)|exp{+z;rf,,} [11] 
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neglecting the effect of the J-coupling during the interval Xi. This assumption is valid in 
the case that the chemical shift frequency difference is much larger than the J-coupling 
(weak coupling approximation). The transformations of the individual state kets are 



i7(A+)|aa> = -|ap> 
l7(A+)|pa> = -i|aa> 
?7(A+)|aP> = -i||3p> 
"-C/(A+)|pp>= |pa> 



which leads to the exchanges of populations sketched in Figure 3c. The population of 
the I ap> state at time point 2 is derived from the thermal equiUbrium population of 
the I acO state. 



The pulse sequence B is given by 90o - X2 - 1 8O90 - T2 - 9O45, with pulse sequence 



delays given by = 



+ /(4J) ^3 = \}(4jj . For the experimental 



results, the delays were set to the values T2 = 44.64ms and T3 = 43.61ms. This sequence 
is designed to separate cleanly the interesting NMR signals deriving from the low-field 
singlet state from the uninteresting signals generated by relaxation of the spin system 
during the second transport interval. 



From standard spin-operator theory, the propagator U(B) for this pulse sequence 
transforms the population operator of state | ap> into a superposition of many operator 
terms, including some that represent antiphase smgle-quantum coherences: 



4 

1 /•_/-._»_ 



[12] 



The (-l)-quantum coherences represented by operators /"/j and iflz generate the 
diagnostic antiphase signals near the chemical shift 62. 
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Sequence B gives rise to no net signals at the chemical shift 82 when applied to in-phase 
Zeeman magnetization, created by Ti relaxation during the second transport interval 
7)^2^. This may be seen fiom the transformations 

C/(5X,C/(5y=-iexpf-^}rr/«--iexpf %]fr/- 

and 

u{BMBf-~i:i; 4-zj/,%- +1/;/- [14] 

Partially-relaxed Zeeman magnetization of spin Ii only generates signals at the chemical 
shift 5i, while partially-relaxed Zeeman magnetization of spin I2 only generates 
multiple-quantum coherences, which do not generate a signal at all. In principle, cross- 
correlated relaxation during the transport interval could also generate spurious signals at 
the shift 82. 

Figure 5 illustrates NMR spectra for a sample of 500 microUtres of a 17.7 mM solution 
of 2,3-dibromothiophene in DMSO-d^ at 20°C. The sample was contained in a 5 mm 
hi^-resolution NMR tube equipped with a Young valve to facilitate degassing. The 
sample contained a water impurity in approximately 27 mM concentration. The sample 
was subjected to three freeze-pump-thaw degassing cycles, each lasting over 30 
minutes, to remove dissolved oxygen. 

The NMR spectra illustrated in Figure 5 were obtained using a Varian Infinity + 
400 MHz NMR system at a field of Bhigh « 9.4 T, using an actively shielded 89 mm bore 
magnet. A standard 5 mm higih resolution NMR probe was used. The 90^ pulse 
duration was around 5 |is. 
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The sample was transported between the high and the low fields using the standaid 
pneumatic sample elevator. The pressure of the elevator was increased during Ihe lift 
operation, and reduced during the insert operation in order to accelerate Ihe transport 
times. Rq)roducible transport times of z£„, = 9±1 jand = 13 ± 2* measured j&om 

when the sample is fiilly loaded into the high field generating coU to when the sample is 
outside the magnet bore. 

The field Biow was estimate using a Hall-effect Gauss meter, with zero field calibrated to 
an accuracy of approximately 10 ^iT using a mu-metal shield far from the magnet. 

For a low or storage field of Blow = 20 mT, the sample was elevated to the top of the 
magnet bore and aUowed to remain there for the time tlf- For even lowra: storage fields, 
the sample was carried physically to calibrated locations within the stray field of the 
magnet, or in the case Blow = 0, a magnetically shielded chamber. 

In the above sample, the proton spin system is defined by 8i= 7.1 1 ppm, 62= 7.72 ppm 
and J = 5.7 Hz. Figure 5a illustrates a conventional proton NMR spectrum, obtained at 
400 MHz in the field of Bugh « 9.4 T. The two doublets of a typical AX spin system 51, 
53 are clearly visible. 

A conventional inversion recovery series (not shown), demonstrates that the spin-lattice 
relaxation constants for the inequivalent sites are very similar in high field (Ti(8i) = 
16.5±0.2s; Ti(82)=17.1±0.2s). 

Figure 5b is a plot of the NMR spectrum for the above sample after it has been 
subjected to the method described with reference to figures 3 and 4. A quasi- 
equiUbrium state was formed by subjected the sanq)le to a hig^i magnetic field and pulse 
sequence A+, then transferring the sauqile adiabatically back to a low storage field. The 
spectrum in figure 3b was taken after the sample was maintained in a low field for a 
time tlf of 1 00s (i.e. much longer than Ti for either 81 or 82). The spectrum of figure 
5b shows diagnostic antiphase signals 55 at around 82=7.72 ppm, and an antiphase 
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dispersion structure 57 around 61 = 7.1 1 ppm. Both of these structures are due to the 
creation of the quasi-equilibrium state in the preparation of the sample. 

The inset of figure 5b shows the region &om. 7.76 to 7.68 expanded. 

The quasi-equilibrium state in the sample is created by manipulating the spins using 
pulse sequence A+ as previously described. To prove that the creation of the singlet 
state is responsible for the features.which exist long after Ti,.a second experiment was 
perfonned. 

The second experiment was performed in the same manner as the first experiment, but 
used the pulse sequence A. instead of the pulse sequence A+. Pulse sequence A. is 
derived &om A+ by changing one of the pulse phases. This exchanges the population of 
state |ap> with that of |pp> instead of with |aa>. The population of |ap> is therefore 
depleted at time point 2 as opposed to enhanced. 

In detail, pulse sequence A. is given by two pulse sequences 90i8o-ti-90,9o. The 
sequence differs fix>m A+ by a 180° phase shift of tiie first pulse. The spin propagator 
for this sequence, under the weak coupling approximation, is: 

U(A_) = exp[-i-|(/^ +/^)jexp{+/;^f^^ } [15] 

which leads to the following ket transformations: 

£/(A.)|aa> = -|aP> 
t;(A.)|pa> = i|pp> 
i7(A.)|ap> = i|aa> 
. £/(A.)|pp>= Ipco 

The population of the | ap> state at tune point 2 is derived fiom the themial equilibrium 
population of the | aa> state. 
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Since the | ap> state is connected adiabatically to the low field singlet, the sequence A. 
as well as A+ can be used to prepare a low field singlet. However, the population of the 
low field singlet deviates fiom the mean in opposite senses for the states formed using 
Af and A.. This variation fi^m the mean in opposite senses means that the diagnostic 
antiphase signals 55 will change sign. 

Thus, in figure 5c, the diagnostic antiphase signals 59 at 62 = 7.72 ppm have a change in 
sign from the equivalent signals 57 in figure 2b. The signals around-5i=7, 1 1 ppm are . 
also perturbed, but in more complicated way since these signals are superpositions of 
several contributions. 

Since the total interval r^^p+r^-^r^^^ = 120^ exceeds Ti by a factor of seven, the 
spectra of figures 5b and 5c prove that the memory of the nuclear spin system has been 
extended significantly beyond the Ti limit 

Figure 6 is a plot of the amplitude of the 7.72 ppm antiphase signals against as fimction 
of tlf for storage magnetic field strengths of OT, 2mT and 20 mT. The amplitude of the 
7.72 ppm antiphase signal is related to the singlet population. 

In a storage field of Biow=1.8 ± 0.5 mT, the time constant for the quasi-equilibriimi 
decay is estimated to be Tqe = 104 ± 5s. The quasi-equilibrium decay time constant had 
no appreciable field dependence for magnetic fields less than Blow = 20 mT. The effect 
of higher magnetic fields on the population of the singlet state will be described with 
reference to figure 5. 

The quasi-equilibrium decay is observed to accelerate when the proton concentration in 
the solution is increased. This indicates that intermolecular proton-proton relaxation is 
a strong contributor to the quasi-equilibrium relaxation, and that highly dilute solutions 
are advantageous to achieve long lifetimes. Quasi-equilibrium relaxation may also be 
caused by paramagnetic impurities such as dissolved oxygen, intermolecular 
interactions with solvent deuterons, and scalar relaxation of the second kmd via the Br 
nuclei. 
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Figures 7a to 7f are plots of simulations of the singlet state population as a function of 
time for a plurality of different fields. La figure 7a Bbw = OT, figure 7b Biow=20mT, 
figure 7c Biow=200mT, figure 5d Biow=500mT, figure 5e Biow=lT and figure 5f Biow=5T. 

The trajectory of the singlet population in a magnetic field is predicted using the 
Liouville-von Neumann equation [4]. For 2,3-dibromothiophene, the distance between 
proton nuclei was estimated to be r = 253 pm. The spinrlattice relaxation in high field 
was assumed to be dominated by intramolecular dipole-dipole relaxation, modelled as 
the isotropic rotational diffusion of a rigid molecule. With this model, the observed Ti 
value of 20s corresponds to a rotational correlation time of Tc=15 ps. 

At low values of Blow, the singlet state population is predicted to be time independent 
over a timescale of hundreds of seconds as seen in figures 7a and 7b, As the storage 
field is increased beyond around 200 mT, a rapid low amplitude oscillation of the 
singlet state population is imposed on a slower decay as can be seen in figures 7c and 
7d. Finally, a higher fields, the amplitude of the oscillations becomes very large and at 
long times the population equilibrates among the four available states, indicating loss of 
all spin order as can be seen in figures 7e and 7f. 

The experimental results show a significant damping of the quasi-equilibrium spin order 
on the timescale of 100s for fields much less than 200mT. This indicates the 
participation of relaxation mechanisms other than intramolecular dipole-dipole 
coupling. Mermolecular dipole-dipole relaxation involving water protons or protons on 
diflFerent DBT molecules is likely to be the major quasi-equilibrium damping 
mechanism. 

In the method described with reference to figures 3 and 4, symmetry operations are 
imposed on or removed fi*om the Hamiltonian by moving the sample fix)m a high to a 
low magnetic field. However, other options are possible. 
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Figure 8 is a flow diagram for performing a method in accordance with an embodiment 
of the present invention where the sample 401 is inherently symmetric in both a high or 
a low magnetic field. An example of such a sample 401 is shown in figure 9. 

A quasi equilibrium nuclear spin ensemble state may be foraied by reacting a sample 
with parahydrogen m a high field environment, the "Pasadena" ejQFect (Russell Bowers 
et al. J. Am Chem. Soc, 109 5541-5542 (1987)), or m a low field envkonment, the 
"Altadena" effect (Pravica et al, Chemical Physics Letters, 145, 255-258 (1988)), 

The symmetric sample 401 is reacted with parahydrogen 403 the presence of a catalyst 
at step S405. This reaction is shown schematically in figure 9. Parahydrogen is a form 
of hydrogen gas in which the hydrogen nuclei are in the nuclear spin smglet state. In 
figure 9, R refers to identical chemical groups. A sample having a quasi equihbiium 
nuclear spin ensemble state is produced 407. 

The sample with the quasi equilibrium nuclear spin ensemble state 407 may then be 
subjected to processes for a time exceeding 3Ti as explained in relation to figure 1. 

To break the symmetry of the sample and thus to allow an NMR signal to be observed, 
the sample is subjected to a photoisomerization process in step S409 by briefly 
illuminating the sample with Ught flash 411. 

A non-equilibrium state of a non-symmetric product 413 is produced as a result of the 
photoisomerization process and the sample is investigated using NMR techniques in 
step S415. 

In the above described method, the sample is kept in a high field environment aU the 
time and symmetry is switched on an off using either a chemical reaction or 
photoisomerization. 
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Figure 10 schematically illustrates a variation on the method described with reference to 
figure 8. Here, the starting sample 501 is an unsaturated non-symmetric molecule 501. 
An example of molecule 501 is shown in figure 11. 

Molecule 501 is reacted with parahydrogen 503 in step S505 in the presence of a 
catalyst and in a low field environment to produce a sample 507 having a quasi 
equilibrium nuclear spin ensemble state. A quasi equilibrium nuclear spin ensemble 
state may be formed by introducmg a sample to. parahydrogen in a low field 
environment, the "Altadena" effect (Pravica et al, Chooadcal Physics Letters, 145, 255 to 
258 (1988)). 

In figure 1 1, A and B refer to non-identical chemical groups. In this case, a quasi 
equilibrium nuclear spin ensemble state is only produced if the external magnetic field 
is low ©aougih that the chemical shift difference between the two hydrogen nuclei is 
negUgible. 

The symmetry is then broken by transferring the sample to a high field to produce non- 
equiUbriiun state 509. An NMR signal may then be detected in step S5 1 1 . 

The flow chart of figure 12 schematically illustrates a variation of the method of figure 
10. In the same manner as figure 10, unsaturated non-symmetric molecule 601 is 
reacted with parahydrogen 603 in step S605 in the presence of a catalyst and in a high 
field environment to produce a sample 607 which is a non-symmetric hydrogenation 
product with enhanced spin order. 

Next, a symmetry operation is irnposed on the Hamiltonian of product 607 by 
transferring the sample adiabatically to a low field environment to form the quasi 
equilibrium nuclear spin ensemble state 609. 

Once the sample 609 has been subjected to the desired processes, symmetry is broken 
by transferring the sample to a high field to produce NMR active sample 611 which can 
be measured in step S613. 
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As previously mentioned, syimnetry operations may be imposed on or removed from 
the Hamiltonian using photoisomerization. 

Chemical reactions, such as photoisomerization, may switch the molecules between 
different conformational or configurational states, in which spin symmetry operations 
are either imposed or removed. As an example, consider the molecule shown in figure 
13a. The two central CH2 groups are chemically, equivalent,, since the molecule has_a 
mirror plane of symmetry. However, the double bonds are readily photoisomerized, 
giving the mixed cis-trans form shown m figure 13b. 

The two central CH2 groups are now chraiically inequivalent, since one has a cis hgand, 
while the other has a trans ligand. A further photoisomerization step generates the all- 
eys form, in which the two central CH2 groups are chemically equivalent again: 

Thus, a Ught flash 701 may therefore be used to switch the symmetry of the spin 
Hamiltonian. 

The flow chart of figure 14 illustrates a fiirther methodin accordance with an 
embodiment of the present invention. As in the case of themethod described with 
reference to figure 3, the sample is not symmetric when in thermal equilibrium in a high 
field environment 801. 

A first KF pulse sequence is applied to the sample in a high field in order to create a 
non-equiUbrium state 803. 

As before, the symmetry of the sample is removed, and due to the populations of non- 
equilibrium state 803, a quasi equilibrium nuclear spin ensemble state 805 is formed 
where the sample 401 is inherentiy symmetric in both a high or a low magnetic field. 
An example of such a sample 401 is shown in figure 9. 
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In this case, the chemical shift differences between the nuclear spins are removed and 
symmetry is achieved by applying a special radiofrequency pulse sequence, a second rf 
pulse sequence, rather than transporting the sample into a low magnetic field. 
Radiofrequency pulse sequences with the appropriate properties are commonly used for 
heteronuclear spin decoupling in solution NMR (see for example, M. H. Levitt, R. 
Freeman and T. A. Frenkiel, Adv. Magn. Reson. 11, 47 (1983)). 

The sample with the quasi equilibrium nuclear spin ensemble state 805 may then be 
subjected to processes for a time exceeding 3Ti as explained in relation to figure 1. 
Throughout this time, symmetry is maintained by the second rf pxilse sequence which is 
continually applied to the sample. The sample is kept in a high field environment. 

To break the symmetry of the sample and thus to allow an NMR signal to be observed, 
the second rf pulse sequence is removed. 

A non-equilibrium state of a non-symmetric product 809 is produced as a result of 
removing the second rf pulse sequence. The sample is investigated using NMR 
techniques in step S811. 
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